Quantum crystal oscillations in a superfluid liquid 
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We study the hydrodynamics of quantum 4 He crystal oscillations in a superfluid liquid with 
involving the dynamics of atomically rough surfaces. It is shown that, due to enhancement of the 
kinetic growth coefficient as the temperature lowers, the reaction force of the liquid applied to the 
4 He crystal changes its character from the inertial to the viscous one vanishing as T — > 0. The 
model is confirmed by the experiments on the oscillations of the 4 He crystal within the temperature 
range 0.54-1.43 K at frequencies 484 and 211 Hz. New type of a hydrodynamic instability is found. 
The instability occurs provided that the oscillation amplitude of the velocity becomes higher than 
~ 3 cm/s. 

PACS numbers: 67.80. -s, 68.45. -v 



INTRODUCTION 

The motion of a helium crystal in the superfluid He 
liquid differs from the motion of the conventional solids 
in their liquid melts. The distinction, in the first turn, 
is associated with the high growth kinetics of the phase 
liquid-crystal 4He interface In fact, the crystalliza- 
tion occurs in the regions with the pressure above the 
phase equilibrium pressure. On the other hand, the melt- 
ing takes place in the regions with the pressure below the 
phase equilibrium pressure. As a result, the shape of a 
crystal and the fluid flow around the crystal can vary 
as the pressure changes beside the crystal. In addition, 
the mass flow across the interface appears and changes 
the boundary conditions and, finally, force applied to the 
crystal. The kinetic growth coefficient K depends on the 
temperature and varies by several orders of the magni- 
tude within the range 0.4 - 1.4 K. This gives a unique 
possibility to study experimentally the effect of phase 
transition dynamics on the hydrodynamics of the motion 
in a superfluid liquid from the case of almost immobile 
crystal surface at high temperatures to the highly mobile 
surface below ~0.55K. Another problem that could be 
solved by setting the crystal into the oscillatory motion 
is associated with the "burst-like growth" effect [H, El- 
The excess pressure at certain temperatures increases the 
growth rate of the crystal facets by several orders of the 
magnitude. This fast-growing state continues to exist 
for several milliseconds after relaxing the pressure to the 
phase equilibrium pressure. It is possible that the peri- 
odic dynamic pressure of the fluid induces such transition 
and stabilizes the anomalous state. 



SURFACE CONTRIBUTIONS TO THE 
DYNAMIC RESPONSE OF A CRYSTAL 

Confining ourselves with the experimental conditions, 
we consider the temperature range 0.4 - 1.42 K, the fre- 
quencies ^200 and ~500 Hz and velocities of the crystal 
motion much less than the speed of sound. In this case 



the superfluid hydrodynamics can be described by two 
independent equations for the superfluid (s) and the nor- 
mal (n) components [4]. In addition, the conditions of 
small oscillations are fulfilled in the experiment. The 
oscillation amplitude A is much smaller than the size 
of the crystal R. This allows us to restrict ourselves 
with the linear terms in the equations of motion. In 
fact, the pressure due to acceleration motion is about 
Space ~ pRa ~ pRAcu 2 , and the Bernoulli pressure is 
SpB ~ pv 2 ~ pA 2 io 2 . For A <C R, the condition 
Space ^ Sps holds for. 

Let us consider the case of a spherical crystal with 
radius R, set into oscillations with frequency u> and rate 
u = to A. The superfluid motion obeys the equations of 
an ideal fluid 

~^ S = V^ S , Aip s = 0, Ps = -p s ^-, (1) 

at 

For the temperatures when p n <c p ~ p s , the situation 
corresponds to the streamline of the sphere with an ideal 
fluid, v ~ v s and p « p s . The boundary condition reads 

it = Tt -P— , V = K6p = K^r-p, Ap = p'-p, 
P PP 

(2) 

where p' is the density of the solid phase, V is the normal 
growth rate of the surface, and K is the growth kinetic 
coefficient. Solving equations Eqs. ((T|) under conditions 
([2]), we obtain a stationary solution for harmonic 1 = 1. 
The force applied to the crystal is given by the expression 
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For high temperatures, coefficient K is small and we ob- 
tain the usual expression for the associated mass. For 
low temperatures and K — > oo, the force of reaction has 
a viscous nature and vanishes. In this limit the move- 
ment of the crystal surface is practically compensated 
by its rcmclting, so that the surface is nearly motionless. 
The harmonics with I > 2 have the following spectrum Q 
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where a is the surface tension which we assume isotropic 
for simplicity. These harmonics cannot be generated by 
the longitudinal oscillations of the crystal. As the tem- 
perature increases, the associated mass of the supcrfluid 
component decreases. The force of reaction in Eq. §3§ 
should be multiplied by a factor p s j p. 

The flow of the normal component along the crystal 
surface contributes to the inertial and to the dissipativc 
parts of the force. The equation for the normal compo- 
nent hydrodynamics reads 



1 dt 



~Vp n + TjAVr, 



(5) 



Following Andreev and Knizhnik |6j, we put the normal 
component velocity to be the same as the interface rate 
V. Thus, the flow of the normal component across the 
boundary is absent, and the crystal can be treated as a 
body with an impermeable surface. The general solution 
of Eq. ([5]) for the oscillations of a sphere in a viscous fluid 
is given in Ref. This expression can be simplified 
since the penetration depth is S vis — y / 2r]/p n uj <C R in 
our experiments. The viscous flow decays in the narrow 
layer beside the crystal interface, while the fluid flow is 
potential in the rest bulk. In this limiting case the force 
acting on the crystal is given by the expression 
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F n = ~YP nR3 ~fa + 3nR 2 y/2r]p n uju. 



(6) 



From Eqs. ^ and ^ it follows that the inertial contri- 
bution depends weakly on the temperature. The decrease 
of the superfluid component contribution is compensated 
by increasing the normal component contribution. A ra- 
tio of the second dissipative term to the first inertial one 
is about ~ 8/R ~ 1CT 3 <C 1. Thus, the total dissipa- 
tive part of the force is determined mainly by the surface 
growth kinetics. 

The above results are applicable to the temperatures 
higher than ~1.2K when the shape of the free growing 
crystal is nearly spherical. The example of such growth is 
given in Ref. [7[. For the further cooling, two roughening 
transitions are observed so that the crystal below ~0.9K 
has a shape of the hexagonal prism. The facets of the 
prism are connected with the atomically rough surfaces. 
The facets have a slow growth kinetics and can be treated 
as immobile in our experiments. Then the crystallization 
and melting will occur at the atomically rough segments 
occupying a small area. For small oscillations, when the 
growth amplitudes of these segments are much smaller 
than their radius of curvature, the estimate of the force 
gives an expression similar to Eq. One should only 
replace the sphere radius with a typical crystal size L and 
the coefficient K should be multiplied by factor £. The 
latter is about a ratio of the area of the active segments to 
the total crystal area. For example, the visual estimate 
of the crystal surface gives us an approximate ratio of 
the rough segments area to the total area of the crystal 
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FIG. 1: Temperature dependence of the oscillator damping. 
The shaded area is a natural damping of the oscillator. The 
explanations are in the text. 



£ ~ 0.1 Since the real shape of the crystal is far from 
spherical, our estimates are rather crude and thus have 
a qualitative character sufficient for our purposes. 

The essential effect of melting and crystallization is 
observed at the temperature for which the condition 
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holds for. For frequencies ~200 and ~500 Hz, the crystal 
size L ~ 0.1cm, £ = 1 and the evaluation with Eq. ([7} 
yields T* = 0.62 - 0.68 K. This is an overestimate since 
the crystal is faceted at such temperature. Provided the 
amount of the atomically rough segments is one tenth of 
the total area (£ = 0.1), the maximum of the damping 
lies within the range T* = 0.49 - 0.54 K, see Fig. [Tj 

For large amplitudes, the regime becomes nonlinear 
since the atomically rough segments are geometrically 
cut by the facets. In the region of the oncoming flow 
the crystallization reduces their area and increases the 
surface curvature. This results in reducing the difference 
between chemical potentials. As a result, the mass flow 
across the interface decreases. On the opposite side of the 
crystal the melting takes place with increasing the cur- 
vature radius of the atomically rough surface segments. 
Since the total volume is closed, the mass flow into the 
liquid results in the total increase of the pressure in the 
container. This pressure should be added to the dynamic 
pressure. The parameters of the influence will depend 
on the size of the crystal, oscillation amplitude, and the 
container volume in which the crystal is grown. The 
crossover to the nonlinear regime occurs at the oscilla- 
tion rates as 
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For T = 0.5 
locity are u* 



0.8 K, the magnitudes of the critical ve- 
= 2 - 200 cm / s. 
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FIG. 2: Left: the first series, T=0.48 K. The crossbar pierces 
the crystal. Right: the second series, T=0.84 K. A crystal 
is located at the platform glued to the loop and (1) tungsten 
tip, (2) a superconducting loop. 



EXPERIMENTAL METHODS 

To set the crystal in motion, we use an oscillator with 
the superconducting U-shaped loop which the crossbar 
is placed in a constant magnetic field ^3 kOe. A tung- 
sten needle is located above the crossbar. The crystal 
nucleates at the needle with applying high voltage which 
induces an additional supersaturation near the surface 
of the tip due to electrostatic pressure. Then the crystal 
remelts down to the crossbar of the oscillator. We use two 
types of oscillators, see Fig. [2] The first represents a su- 
perconducting wire with the resonance frequency 484 Hz. 
The second is made with a glassy table glued with epoxy 
to the crossbar in order to reduce the frequency to 211 Hz. 
The first design is suitable for the measurements below 
0.9 K, when the crystal is cut by the facets and remains 
at the crossbar for a long time. The second design of the 
oscillator (see Fig. 2, right) prevents from the remelting of 
the crystal in the hydrostatic pressure gradient and from 
descending the crystal to the bottom of the container. 

The latter design can be used within the whole temper- 
ature range. The oscillations of the loop are generated 
with an alternating current. The voltage arising due to 
the motion of a crossbar in the magnetic field is amplified 
and registered with the computer-supplied measurement 
system. We use two methods. In the first one the am- 
plified voltage of the loop is fed to the digital amplitude 
stabilizer Q and transferred after the phase shifter and 
current amplifier back to the loop. The self-exciting oscil- 
lator allows us to measure continuously the frequency and 
the damping of the system at the fixed oscillation ampli- 
tude. The absolute magnitudes of the decrement and fre- 
quency are determined from the free damped oscillations 
after the break-up of feedback. The second method repre- 
sents a measurement of the amplitude-frequency charac- 
teristics of the system. This method is employed for large 
crystals yielding high attenuation and frequency shift in 
the system. 

In the experiment we use the oscillation mode corre- 



sponding to the synchronous oscillations of the vertical 
parts of the loop. The horizontal crossbar moves as a 
whole, producing an additional inertial force at the free 
ends of the vertical segments of loop. The equation de- 
scribing this mode is derived from the equation for the 
vibrations of the rod with one fixed end [9]. For deriva- 
tion, we augment the inertial force resulted from the finite 
mass of a half of the crossbar, exerted to the free end of 
the vertical segment of the loop. Finally, we arrive at 

cos(fcL) coth(fci) + 1+ 
+akL[cos(kL) sinh(fcL) — s'm(kL) coth(fcL)] = 

r z E 2 



PP 



a = -^—(Mi + M c + M c — - ) 

M .i v 2p'l+icoT J 



(9) 



Here L is the length of the vertical segments of the loop, r 
is the radius of the wire, E is its Young modulus, M 0j i are 
the masses of the vertical segments and horizontal cross- 
bar, M c is the helium crystal mass, and R is the radius 
of the crystal. An additional bending stress at the verti- 
cal ends, associated with the rotational oscillations of the 
crossbar, changes slightly the main mode frequency. In 
general, there exist other oscillation modes, for example, 
the torsional mode with the out-of-phase motion of the 
vertical parts of the loop or the modes with exciting the 
proper oscillations of the crossbar. These modes have 
higher frequencies as compared with the principal mode 
described by Eq. ([9]) . In our experiment these additional 
modes were not excited. 

The damping rate is determined by the imaginary part 
of the frequency u) found from solving Eq. (J9j) . We define 
the decrement of the system S as the logarithm of a ratio 
of two oscillation amplitudes, one at the time moment t 
and the second at the time t + T later by the period T , 
i.e., 5 = 27rIm(u;)/Re(u;). 

Using the data on the temperature behavior of the ki- 
netic growth coefficient, in Fig. Q] we show the calcula- 
tion for the attenuation of the composite oscillator with a 
crystal of radius 1 mm. Above the first roughening tran- 
sition temperature ^1.25 K the decrement of the system 
at the both frequencies is less than 10~ 3 . This value is 
smaller than the damping of the free loop. Below the 
both roughening transition temperatures we use Eq. Q 
with a correction of the growth kinetic coefficient by a 
factor of £ = 0.1, see The solid line is a calculation 
for the first series of experiments and the dashed line is 
for the second experimental series. The increase of the 
crystal size shifts the dependence S(T) towards higher 
temperatures. 

The curves Fast A, B are calculated by using the crys- 
tal growth rates in the abnormal state [l3| for the fre- 
quencies of the first and second runs, respectively. The 
crystals have the same size and mass but dependencies 
S(T) at Fig. [1] differ by the oscillation frequency. From 
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the plot one can see that the transition of a crystal to 
the "burst-like growth" regime, which is accomplished 
by the drastic acceleration of the growth rate for all crys- 
tal facets, results in a jump of decrement. The observa- 
tion of such jump with keeping new magnitude of damp- 
ing would evidence for the transition of the crystal to a 
steady anomalous state. 



EXPERIMENTAL RESULTS 

The frequency and damping of the free oscillator in su- 
perfluid helium at low oscillation amplitudes are weakly 
dependent on the temperature. The frequency increases 
by 0.01% along the melting curve for cooling from 1.8K 
to 0.48 K. The decrement is governed by the sound emis- 
sion into the mounting base and lies within the range of 
0.002 - 0.003. As the oscillation amplitude increases to 
the velocity ~1 m/s, the damping enhances by one or- 
der of the magnitude. For the small amplitudes, when 
the maximum rate does not exceed ~10 cm/s, the decre- 
ment increases approximately linearly. 

The experiments are performed at 1.24 K and 1.43 K as 
well as below the first and second roughening transitions 
(0.54 K, 0.64 K, 0.74 K and 0.84 K). Placing a crystal at 
the crossbar reduces the free oscillation frequency of the 
oscillator in accordance with the model. An additional 
damping contributed by the crystal above the roughen- 
ing transition temperatures does not exceed the natural 
damping of the oscillator in the agreement with the cal- 
culation after Eq. (j9|) . The qualitative behavior of damp- 
ing below the roughening transitions agrees well with the 
model proposed. In Fig. [3] the results are shown for the 
measurements of damping at the crystal before, during 
and after growth at 0.733 K. Before starting the growth 
the crystal of diameter 2.4 mm and thickness 0.65 mm 
has the large rounded edges giving an additional signif- 
icant damping. With the beginning of the growth the 
pressure in the fluid increases, entailing the increase of 
the edge curvature and decrease of their area. As a result, 
the damping of the system drops to the magnitude close 
to the damping of the free loop. As the crystal growth 
stops, the pressure in the fluid relaxes to the phase equi- 
librium pressure, the size of the edges increases, and the 
attenuation intensifies. In Fig. [1] this evolution is shown 
by the squares in the small amplitude limit. The starting 
point is close to the theoretical curve for the first series 
of experiments. The drop of damping during the growth 
is shown by the downward arrow. The recovery to the 
higher magnitudes of the decrement, shown by the up- 
ward arrow, agrees also with the model proposed. The 
reason is that the increase of the crystal size shifts the 
dependence S(T) to the high temperature region. On 
the right-hand branch of the curve, this leads to the in- 
creased attenuation, as is observed in the experiment. 
The normal frequency of the system falls gradually from 
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FIG. 3: The size effect of the atomically rough edges on the 
oscillator damping. The damping in the system: (A) before 
growth, (•) during growth, (□) at the end of growth. The 
solid line is the free loop attenuation 



442 Hz to 384 Hz with accumulating the crystal mass in 
the growth process. 

The enhancement of the oscillation amplitude leads to 
the appearance of instability which results in sliding the 
crystal down from the oscillator crossbar. The instability 
is observed in the both series of experiments at a max- 
imum oscillating rate more than ^3 cm/s. For the first 
series of experiments in which the crystal is pierced by 
the crossbar, the crystal starts to grow intensively until 
the crystal drops to the bottom of the cell. In the sec- 
ond run the crystal glides drastically down as the criti- 
cal rate is reached. Expecting that the crystal could be 
attracted to the oscillator as it is observed for the oscil- 
lations of the crossbar beside the surface (Fig. |4}. We 
have attempted to put the crystal from the needle to the 
crossbar oscillating with the large amplitude. However, 
that produced no result. For the oscillating rates with a 
maximum velocity above ^3 cm/s, the crystal bounces 
from the crossbar. The existence of such instability lim- 
its the maximum magnitude of the overpressure by ^0.2 
mbar at the critical point of the liquid streamline of the 
crystal. This value is one order of the magnitude less 
than that required for the transition of the crystal into 



the abnormal state 11 



Below roughening transition temperatures Tri^ in 
the steady state the crystal is suspended at the needle 
tip for hours jzj. The Rayleigh- Taylor instability does 
not develop at the bottom facet during this time as a re- 
sult of a giant potential barrier which should be overcome 
via either thermal activation [12j or quantum tunneling 
[l3j |. For this reason, the instability observed in our ex- 
periment is associated directly with the motion of the 
crystal. 

In the reference frame associated with the oscillator 
crossbar at the acceleration maximum one half of the 
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FIG. 4: The attachment of the atomically rough crystal sur- 
face to the oscillating crossbar. The temperature is 1.39 K, 
crossbar rate swing is 20 cm/s, and the spacing between the 
crossbar and undisturbed crystal surface is 0.7 mm. 



as well. The model proposed describes well the quali- 
tative features of the phenomenon and satisfactory the 
numerical magnitudes of the damping. Unfortunately, 
the appearance of hydrodynamic instability for the oscil- 
lations at a rate higher than 3 cm/s prevents us from ob- 
taining high overpressures at the crystal surface. Hence, 
the questions whether it is possible to prepare the helium 
crystal in the stationary abnormal state by oscillating the 
crystal has remained open. 
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crystal is under conditions of appearing the Rayleigh- 
Taylor instability. However, the acceleration magnitude 
is about ~ 10 4 cm/s 2 ~ lOg at the velocity ~ 3 cm/s and 
frequency ~ 500Hz, g being the acceleration of gravity. 
This is still insufficient for overcoming the potential bar- 
rier during the oscillation half-period [H, . An accel- 
erating motion can give rise to the Richtmayer-Meshkov 
instability [l4[. However, the estimates show that the 
acceleration magnitude is too small to develop this in- 
stability as well. 

The fluid flow along the facet, as is shown by Andreev 
15{ , decreases the energy of the steps normal to the fluid 
flow and increases the concentration of steps. These pro- 
cesses, in principle, may favor the classical mechanisms 
of the facet growth. 

To conclude, the origin of the instability observed as 
sliding the crystal from the oscillator crossbar has not 
been established. 



SUMMARY 



The kinetics of atomically rough surfaces in a faceted 
4 He crystal affects significantly the motion of the crystal 
in a superfluid liquid and governs the hydrodynamics of 
the streamline. The damping of the oscillating crystal is 
determined mainly by the temperature behavior of the ki- 
netic growth coefficient and the total area of rough edges 



* vlt@isssph.kiae.ru 
[1] H. Alles, S. Balibar and A. Ya. Parshin, Rev. Mod. Phys. 

77, 317 (2005). 
[2] V. L. Tsymbalenko, Phys. Lett. A 211, 177 (1996). 
[3] P.J. Hakonen, J.S. Penttila, J. P. Ruutu, J. P. Saramaki, 
G. Tvalashvili, A.V. Babkin and A.Ya. Parshin, Phys. 
Rev. Lett. 76, 4187 (1996). 
[4] L. D. Landau and E. M. Lifshits, Fluid Mechanics, (Perg- 

amon Press, 1987). 
[5] S. N. Burmistrov, L. B. Dubovskii and V. L. Tsym- 
balenko, Phys. Rev. E 79, 051606 (2009) 
[6] A. F. Andreev, V. G. Knizhnik, Zh. Eksp. Teor. Fiz. 83, 

416 (1982) [Sov. Phys. JETP 56, 226 (1982)]. 
[7] V. L. Tsymbalenko, Fiz. Nizk. Temp. 21, 162 (1995) [Low 

Temp. Phys. 21, 120 (1995)]. 
[8] V. L. Tsymbalenko, IET 3, 129 (2009). 
[9] L. D. Landau and E. M. Lifshits, Theory of Elasticity, 
(Pergamon Press, 1968). 
[10] V. L. Tsymbalenko, Zh. Eksp. Teor. Fiz. 119, 1182 (2001) 

[JETP 92, 1024 (2001)]. 
[11] V. L. Tsymbalenko, J. Low Temp. Phys. 121, 53 (2000). 
[12] S. N. Burmistrov, L. B. Dubovskii and V. L. Tsym- 
balenko, J. Low Temp. Phys. 162, 383 (2011). 
[13] S. N. Burmistrov, L. B. Dubovskii, J. Low Temp. Phys. 

150, 295 (2008). 
[14] S. N. Burmistrov, L. B. Dubovskii and V. L. Tsym- 
balenko, J. Low Temp. Phys. 162, 391 (2011). 
[15] A. F. Andreev, Zh. Eksp. Teor. Fiz. 106, 1219 (1994) 
[Sov. Phys. JETP 796, 660 (1994)]. 



